EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH 



CERN-EP/2000-145 
November 24, 2000 



O 
O 
O 

(N 

O 
D 

Q 



Search for Neutral Higgs Bosons of the Minimal 
Supersymmetric Standard Model in e + e _ Interactions at 

= 192 - 202 GeV 



The L3 Collaboration 



> 

o 

(N 

o 
o 

X ' Abstract 

<D ' 

A search for the lightest neutral CP-even and the neutral CP-odd Higgs bosons 
of the Minimal Supersymmetric Standard Model is performed using 233.2 pb _1 
of integrated luminosity collected with the L3 detector at LEP at centre-of-mass 
energies 192 — 202 GeV. No signal is observed and lower mass limits are given as a 
function of tan/3 for two scalar top mixing hypotheses. For tan/3 greater than 0.8, 
they are > 83.4 GeV and ttt-a > 83.8 GeV at 95% confidence level. 
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1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) Tj requires two Higgs doublets. This 
gives rise to five Higgs bosons: a charged scalar pair, two neutral CP-even, the lightest of which 
is called h, and a neutral CP-odd, A. The two most important production mechanisms in e + e~ 
collisions are: 

e + e"^Z*^hZ (1) 
e+e--^Z*^hA. (2) 

The cross section of the process (0) is smaller than that for the similar production of the 
Higgs boson in the Standard Model. This process is dominant at low values of tan/3 (tan (5 < 5), 
where tan/3 is the ratio of the two Higgs vacuum expectation values. The pair-production of 
Higgs bosons (j2J) takes over at high values of tan/3. 

Previous searches for the h and A bosons were reported by L3 |2j and other experiments j^j. 
In this paper, we present the results of the search for the h and A bosons using the data collected 
with the L3 detector jlj at centre-of-mass energies yfs = 191.6 — 201.7 GeV, corresponding to 
233.2 pb~ of total integrated luminosity. The sensitivity to the production of MSSM neutral 
Higgs bosons is improved by combining the results of these analyses with our previous searches. 

2 Data and Monte Carlo Samples 

In 1999 the L3 detector collected data at LEP at average centre-of-mass energies ^/s = 
191.6 GeV, 195.5 GeV, 199.5 GeV and 201.7 GeV, corresponding to the integrated luminosi- 
ties 29.7 pb~\ 83.7 pb~\ 82.8 pb" 1 and 37.0 pb^ 1 respectively. 

The cross sections of processes Q and (J2J) and the decay branching ratios of h and A are 
calculated using the HZHA generator [H] . For efficiency studies, Monte Carlo samples of Higgs 
events are generated using PYTHIA |H| and HZHA. 2000 Monte Carlo events are generated for 
each mass hypothesis. For hA samples the masses, and ??2a, of the h and A bosons range 
from 50 to 95 GeV in steps of 5 GeV. For hZ samples is chosen in steps of 5 GeV from 50 to 
95 GeV and in steps of 1 GeV from 95 to 110 GeV. For the background studies, the following 
Monte Carlo programs are used: PYTHIA (e + e~ — > qq(7), e + e~ — > TIL and e + e~ — > Ze + e~), 
KORALW [Zj (e + e~ -> W+W"), KORALZ [B] (e+e" -> r+r"). Hadron production in two- 
photon interactions is simulated with PYTHIA and PHOJET 0. EXCALIBUR PU| is used 
for other four fermion final states. The number of simulated background events for the most 
important background channels is more than 100 times the corresponding number of expected 
events. 

The L3 detector response is simulated using the GEANT program [IT] , which models the 
effects of energy loss, multiple scattering and showering in the detector. The GHEISHA pro- 
gram ^2] is used to simulate hadronic interactions in the detector. Time dependent inefficiencies 
are also taken into account. 

3 Event Selection 

For the hA production, the following decay modes are considered: hA — > bbbb, hA — > bbr + T~ 
and hA — > r + r~bb. In the case of hZ, four event topologies covering approximately 98% of 
possible final states, are considered: qqqq, qqz/z/, qql + l~(l = e,fi,r) and T + r~qq. The searches 
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in channels with hadronic decays of the h boson are optimised for the dominant h — * bb decay 
channel. The analyses qquu and qql + l~(l = e, fi) are the same as devised for the Standard 
Model Higgs search [THj . 

A common selection is applied to both the hA and hZ searches. In the four-jet channel, it 
mainly reduces the two-photon interaction background while keeping the signal efficiency high, 
then a neural network is used to build a discriminating variable. In the tau channel, first a 
selection is devised, then an optimal variable based on a likelihood approach is defined. 

3.1 hA — ►bbbb and hZ — ► bbqq Selection 

The signature of both the hA— >bbbb and hZ— >bbqq final states is four hadronic jets and the 
presence of b-hadrons. The dominant backgrounds come from qq(7) production and hadronic 
decays of W and Z pairs. 

High multiplicity events are selected and their visible energy, E vis , is required to be greater 
than 0.6y/s and less that lA^/s. Events with perpendicular imbalance greater than 0.35E vis 
or with a lepton whose energy exceeds 65 GeV are rejected to suppress semileptonic W pair 
decays. Initial state radiation events are further suppressed by requiring P^ is / (m vis -mz) < 0.4, 
where P^ is is the longitudinal component of the missing momentum, m vls is the visible mass 
and mi the Z boson mass. The remaining events are then forced into four jets using the 
DURHAM algorithm ^3] and a kinematic fit requiring energy and momentum conservation 
(4C) is performed. The number of expected and observed events in the data, together with the 
signal efficiencies for the selection cuts are listed in Table 

After the selection, discriminating variables are combined into a feed forward neural net- 
work with one hidden layer and three output nodes. The inputs include the probability 
that the jets contain b-quarks ^H], hereafter termed Bxag, event shape variables and mass in- 
formation. The information on the event shape includes the event sphericity, the value of the 
DURHAM jet resolution parameter for which the event is resolved from three to four jets, the 
longitudinal component of the missing momentum and the event thrust. The mass information 
is summarised in a x 2 variable defined as: 

X 2 (m x , m h ) = (£|f - (m x + m h )) V + (A f n - \m x - m h \) 2 w A , (3) 

where m x is either mz or m&. The pairing used gives the minimum value for the difference 
squared (Ay — \m x —m h \) 2 between the measured and expected dijet mass differences; E™ 111 is 
the corresponding sum of the dijet masses. The weights, and w&, are derived from the mass 
resolutions, and their ratio is 3/5. In this way, the shape of the neural network output is made 
almost independent of the mass hypothesis. The polar angle of the Higgs boson, 9, is also used 
as input for the hA analysis. It gives additional separation between the hA signal and W + W _ 
background, due to the different spins of the W and the Higgs bosons. The distributions of the 
discriminating input variables are shown in Figure H 

The Neural Network has three output variables which correspond to the signal, O Higgs, the 
qq(7) final state, O qq , and the W + W~ final state hypotheses, Oww- The discriminating variable 
is then obtained from the combination NN = O Higgs x (1 — O qq ) x (1 — Oww)- Two different 
neural networks are used for the hA and hZ analyses, and the events are classified as hA or hZ 
according to the largest value of the discriminating variable. The neural network for the hZ is 
optimised for m h = 100 GeV at ^/s = 192-196 GeV and m h = 105 GeV at y/s = 200-202 GeV. 
For the hA search, the neural network is optimised for = mx = 85 GeV. Examples of the 
discriminating variable distributions of the two analyses are shown in Figure El Good agreement 
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is observed between the data and the expected background. An independent analysis based on 
a likelihood approach [T^j validates the present results. 

The highest signal sensitivity region corresponds to large values of NN. For illustrative 
purposes, in Table we report the number of the observed and expected events selected after 
the cuts NN >0.9 for the hA analysis, and NN >0.5 for the hZ analysis. The mass distributions 
for the hA search after the cut NN >0.9 is shown in Figure EH- 





192 GeV 


196 GeV 


Cut 


Selection 


NN hA >0.9 


iViV h z>0.5 


Selection 


NN hA >0.9 


NN hZ >0.5 


Data 


430 


2 


4 


1186 


3 


33 


MC 


433.6 


1.0 


9.5 


1197.5 


3.2 


26.0 


qq 


201.2 


0.4 


3.0 


531.3 


1.2 


7.1 


W + W~ 


217.9 


0.2 


4.2 


622.3 


0.7 


11.9 


zz 


14.5 


0.4 


2.3 


43.9 


1.3 


7.0 




95.0% 


42.0% 


65.9% 


95.0% 


42.0% 


65.9% 


£hZ 


93.7% 


12.6% 


52.7% 


93.7% 


12.6% 


52.7% 




200 GeV 


202 GeV 


Cut 


Selection 


NN hA >0.9 


NN hZ >0.5 


Selection 


NN hA >0.9 


NN hZ >0.5 


Data 


1198 


2 


9 


506 


1 


6 


MC 


1118.4 


2.0 


9.1 


507.8 


1.0 


4.8 


qq 


478.3 


0.7 


2.5 


215.0 


0.3 


1.5 


W+W~ 


595.4 


0.3 


3.3 


271.8 


0.2 


1.7 


zz 


44.7 


1.0 


3.3 


21.0 


0.5 


1.6 


£hA 


94.1% 


37.9% 


37.9% 


95.9% 


38.3% 


40.4% 


£hZ 


91.0% 


5.7% 


40.9% 


92.3% 


5.0% 


36.3% 



Table 1: Number of events observed and expected in the four-jet channels, after the selection 
and after cuts on the discriminating variables, NN^a and NN^ Z . The hA signal efficiencies are 
quoted for m\ = = 85 GeV. The hZ signal efficiencies correspond to = 100 GeV at 
= 192 - 196 GeV and m h = 105 GeV at = 200 - 202 GeV. 



3.2 hA— >bbr + T , hZ— >bbr + T and hZ — >t + t qq Selections 

The signatures of hA— >hhr + T~ 0, hZ^bbr + r _ or hZ^r + T~qq events are a pair of taus 
accompanied by two hadronic jets. For each of the channels hA and hZ an analysis is optimised 
based either on the tau identification or on the event topology by requiring four jets with two of 
them being narrow and of low multiplicity. The main background results from W-pair decays 
containing taus. 

The hZ analysis is similar to the one described in detail in Reference [Tj3 The hA selection 
is optimised for lower Higgs masses by omitting the cuts on the opening angles of the jets and 
tau pairs and on the invariant mass of the tau pair, m TT . The invariant mass of the hadronic 
jets, m qq , must be between 5 GeV and 125 GeV. The ratio of the sum of the energies of the 
tau decay products over the sum of the energies of the jets is required to be less than one and 
the value of the missing momentum vector in the rest frame of the Higgs must be less than 

^Both the decay modes (h— * bb, t + t~) and (h— > t + t~, bb) are considered. 
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40 GeV. Finally, the cosine of the polar angle of the Higgs boson, | cosG|, has to be less than 
0.8. The number of events observed, the number expected from background processes, and the 
signal efficiency for the two selections are listed in Table El 

For each event class j (ZZ, W + W~, qq, Ze + e~, hA, hZ), a probability function, /j, is 
constructed, where i denotes the variables considered. These are the Bxag for each hadronic 
jet, m qq and m TT . They are presented in Figure |H The probability, p^, of an event to belong 
to class j, based on the value of the variable i, is then defined as 



Pi 



n 



(4) 



Finally, the probabilities for the individual variables are combined by calculating the likelihood 
that the event belongs to the either signal class: 



FhA 



J2k YliPk 



and 



J2k Yli p\ 



(5) 



Events retained by both the hA and the hZ selections, are classified according to highest 
value of the likelihood, as hA or hZ candidates. An example of the distribution of the discrim- 
inating variable for the hA search is shown in Figure Good agreement between the observed 
data and the expected background is found. The mass distribution for the hA search after an 
additional cut on the Bxag is shown in Figure Eb- 
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Table 2: Number of events observed and expected in the tau selection. Efficiencies for the hA 
signal are quoted for = mj, = 85 GeV. For the hZ signal, they are quoted for m\ l = 100 GeV 
at y/i = 192 - 196 GeV and for m h = 105 GeV at y/s = 200 - 202 GeV. 



4 Results and Interpretation 

A good agreement between data and expected background, both in the total number of events 
and in the shape of the distributions, is observed in all the channels analysed. The mass 
distributions of the events in the highest sensitivity region are not compatible with a signal for 
any mass hypothesis. Therefore no evidence of the production of the h and A bosons is found. 
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The results of the search for the hA and the hZ production are interpreted in the framework 
of the constrained MSSM (CMSSM) assuming unification of the scalar fermion masses, unifica- 
tion of the gaugino masses and unification of the trilinear Higgs-fermion couplings at the GUT 
scale. This choice has little impact on the phenomenology of the Higgs bosons but reduces 
significantly the number of free parameters. The remaining free parameters are tan/3, m A , the 
gaugino mass parameter, M2, the universal scalar fermion mass, mo, the common scalar quark 
trilinear coupling, A, and the Higgs mixing parameter, //. 

Two benchmark scenarios jTHj are considered. In the first one, termed "maximal mixing" , 
the CMSSM parameters are chosen such that acquires its maximal value for any given value 
of m A and tan/3. The second scenario corresponds to vanishing mixing in the scalar top sector 
and is referred to as "minimal mixing" . 

The CMSSM parameters are chosen as follows: uiq = 1 TeV, \x = —200 GeV, M2 = 
200 GeV. The mass of the top quark is fixed to 175 GeV. The maximal mixing scenario is 
realised at X t = A — fj, cot (3 = TeV, where X t is the parameter which controls the mixing 
in the scalar top sector. The minimal mixing corresponds to X ( = 0. Keeping these values 
fixed, a scan over the two remaining independent parameters, tan/? and toa, is performed in 
each mixing scheme in the ranges: 0.5 < tan/5 < 30 and 10 GeV < tua < 1 TeV. 

To set exclusion limits on the CMSSM parameters, the confidence level, CL, that the ex- 
pected signal is absent in the data, is calculated ;20j for each point (tan/3, m A ) of the scan. The 
full distributions of the discriminating variables, NN and F, are used in this calculation. 

Systematic and statistical uncertainties on the signal and background are evaluated using 
the same procedure as in the Standard Model Higgs search The main sources of systematic 
uncertainties are detector resolution, selection procedures, theoretical uncertainties and Monte 
Carlo statistics. The overall systematic uncertainties are estimated to be 4% on the predictions 
for the expected signal and 10% for the background events. Bins of the final variables with 
a signal-over-background ratio in the Monte Carlo of less than 0.05 are not considered in the 
calculation of the CL. This cut is chosen to minimise the effect of systematic uncertainties on 
the average CL as calculated from a large set of Monte Carlo experiments. 

The results of the MSSM Higgs search at lower y/s [2] are combined with those presented 
in this paper. Figure El shows the region of the (tan/3, m^) plane and (tan/3, m A ) plane excluded 
by L3 for the maximal mixing and minimal mixing scenarios. 

For the CMSSM parameters considered and assuming tan/3 greater than 0.8, this results in 
lower mass limits at the 95% CL of: 

m h > 83.4 GeV, m A > 83.8 GeV, 

which compare to the median expected limits in the absence of a signal of > 85.6 GeV and 
m A > 85.7 GeV. 

The exclusion plots for the minimal mixing scenario present a small unexcluded area in the 
low tan/? region at low values of m A where the decay h — > AA is allowed but is not investigated 
among the signatures described above. 

For 0.8 < tan/? < 1.8 values of m A up to 1 TeV are ruled out for any mixing scenario 
allowing to exclude this tan/3 region in the CMSSM, for the top mass < 175 GeV. 
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Figure 1: Distributions of a) Bxa g , b) logarithm of the \ 2 variable, c) cosine of the polar angle 6 
of the Higgs boson and d) sphericity after the four-jet selection. The points represent the data 
at \fs = 192 — 202 GeV, the open histograms are the expected Standard Model backgrounds, 
and the hatched histograms are the expected hA — > bbbb signal scaled by a factor 50, for 
= ttia = 85 GeV at tan/5 = 30. 
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Figure 2: Distributions of the discriminating variable in the four jet channel for a) hA — > bbbb 
and b) hZ — > bbqq. The points show the data collected at ^/s = 192 — 202 GeV, the open 
histograms are the expected Standard Model backgrounds and the hatched histograms are the 
expected signals scaled by a factor 10. The discriminating variables are constructed assuming 
equal mass hypothesis m h = m A = 85 GeV at tan/5 = 30 in a) and m h = 100 GeV at tan/3 = 1 
in b). 
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Figure 3: Distributions of a) the sum of the reconstructed dijet masses in the hA^bbbb 
channel after the cut NN >0.9, b) the sum of the dijet and the ditau masses in the hA— >bbr + T~ 
channel after a cut on the Bxag- The points represent the data collected at ^/s = 192 — 202 GeV, 
the open histogram is the expected Standard Model background and the hatched histogram is 
the expected signal for m h = = 85 GeV at tan/5 = 30. 
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Figure 4: The distributions for the hA— >bbr + r~ search of a) the BT ag for the most energetic 
hadronic jet and b) the least energetic hadronic jet, c) the reconstructed mass for the hadronic 
system, and d) the reconstructed mass for the leptonic system. The points are the y/s = 
192 — 202 GeV data, the open histograms the expected Standard Model background and the 
hatched histograms the expected hA— >bbr + r~ signal for m h = = 85 GeV at tan/? = 30. 
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Figure 5: Distribution of the discriminating variable of the hA^bbr + r _ selection in the 
hypothesis = rriA = 85 GeV at tan/5 = 30. The points are the data collected at 
yfs = 192 — 202 GeV, the open histogram is the expected Standard Model background and 
the hatched histogram is the expected signal for the hA search multiplied by a factor 10. 
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Figure 6: Exclusion plots in the (tan/3,mh) and (tan/?,mA) planes at the 95% CL for the minimal 
and maximal mixing scenarios. The hatched area represents the exclusion and the crossed area 
is not allowed by the theory. The horizontal hatched area corresponds to < 10 GeV and 
was previously excluded at LEP 
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